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SUMMARY 


An analytical study has been performed to investigate excessive heating in the 
tile-to-tile gaps of the Shuttle orbiter thermal protection system (TPS) . The exces- 
sive heating caused visible discoloration and charring of the filler bar and strain 
isolator pad (SIP) used in the attachment of tiles to the aluminum substrate. The 
study concentrated on determining if tile-to-tile steps and gaps measured after 
flight could create a sufficient pressure disturbance around a stepped tile to cause 
hot gas flow in the gaps during entry. Techniques based on existing experimental 
data were developed to predict the local pressure disturbances. In addition, a tech- 
nique was devised to approximate the energy potential of the air flowing from the 
local boundary layer into the tile-to-tile gap. A tile internal flow model, devel- 
oped previously for the pressure loading on a tile during ascent, was modified to 
account for compressibility in the flow in the tile-to-tile gaps during entry. The 
flow model predicts mass flow rates and pressure distributions in the tile-to-tile 
gaps, in the SIP, and within the tiles. It utilizes analogies between conventional 
heat conduction and internal gas flow in the calculation procedure. With the flow 
rates established, a thermal model of adjacent tiles, the tile-to-tile gap, filler 
bar, and SIP is used to predict temperature distributions throughout the system. An 
energy balance is maintained on the air flowing through the tile system. 

These analytical tools are used to predict conditions under which damage to the 
filler bar and SIP could occur. 


INTRODUCTION 

The Space Transportation System (STS) orbiter employs silica tiles in the design 
of a reusable thermal protection system (TPS). (See fig. 1(a).) The tiles are 
typically 6 in. by 6 in. and are of variable depth, depending on their location on 
the orbiter. The tiles are attached to the aluminum substrate with a strain isolator 
pad (SIP) made of high-temperature nylon fibers that is bonded to the bottom of the 
tiles and the aluminum with silicone rubber. The bonded area of the SIP is 5 in. by 
5 in. on the bottom of the tiles- The aluminum substrate outside the tile/SIP/ sub- 
strate bonding area and directly under the tile-to-tile gaps is protected with a 
high-temperature nylon fiber filler bar that has been heat-treated to withstand 
higher temperature exposure than SIP. The filler bar is bonded to the aluminum sub- 
strate, but it is not bonded to the bottom of the tiles. The filler bar also has a 
silicone rubber membrane that is applied to the top of the bar to restrict gas flow 
under the tiles. During the first Shuttle flights, the filler bar experienced higher 
than expected heating at random locations on the lower side of the fuselage and wing. 
The excessive heating on the filler bars caused varying degrees of damage - from 
slight discoloration of the silicone rubber membrane to complete charring of the 
membrane and high-temperature nylon fiber. The damage to the filler bars occurred on 
the bottom of the fuselage and wings in geometrically flat regions where the aero- 
dynamic pressure gradient is essentially zero. Ceramic fiber gap fillers are used 
between tiles in regions where pressure gradients exist, such as near the leading 
edge of the wings. Gap fillers were not used at these damaged filler-bar locations. 

Postflight inspection of the vehicle and its TPS yielded no obvious rationale 
for the char pattern exhibited in figure 1(b). All tiles were intact and firmly 



attached. The tile step and gap measurements were somewhat irregular and often 
exceeded preflight measurements. The mismatches were as large as 0.120 in. in step 
height and 0.160 in. in gap width. Attention was focused on the tile-to-tile steps 
as a possible cause of the filler-bar heating, since local pressure disturbances 
would be created on the upstream and downstream faces of the stepped tile. A com- 
bination of high and low pressures around a tile causes airflow in the tile-to-tile 
gap, in the SIP, in the filler bar, and in the tile itself. This airflow, at a high 
energy potential during entry, can cause damage to the filler bar and SIP if they are 
exposed to these airflow rates over a sufficient period of time. It was not immedi- 
ately obvious that the observed tile step heights were large enough to cause signifi- 
cant airflow in the tile-to-tile gaps, since they were deeply submerged in the local 
boundary layer during entry. Boundary-layer thicknesses on the forward part of the 
wing and fuselage were approximately 2.8 in. and 4 in. , respectively, during peak 
entry heating. An analytical effort was undertaken to define the step heights and 
gap widths which could produce local conditions capable of charring the filler bars 
to the extent observed. 

This report describes the techniques used and the results of an analytical 
investigation of excessive heating in the gaps between high-temperature reusable 
surface insulation (HRSI) tiles of the thermal protection system (TPS) on the lower 
surface of the fuselage and wings of the Shuttle orbiter. The investigation centered 
around determining analytically whether the observed postflight tile-to-tile steps 
and gaps were large enough to generate sufficient hot-gas flows in the gaps to cause 
damage to the filler bar and SIP. Techniques were developed for estimating the pres- 
sure disturbances due to a stepped tile and for calculating the disturbance— induced 
mass flow rates in the tile-to-tile gaps, the filler bar, the SIP, and the tile 
itself. These results were used to perform a thermal analysis in the tile-to-tile 
gap to determine temperature response to the hot-gas flow. Calculations were per- 
formed at one location on the fuselage and one location on the left wing, where dam- 
aged filler bars were observed on the first flight. The results of the analysis are 
presented in this report. 


SYMBOLS 


A cross-sectional area for flow 

c 

A area for pressure loading of SIP 

Pl 

A area for pressure loading of filler bar 

P2 

Ag surface area 

A area of SIP which is stressed 

A area of filler bar which is stressed 

^2 

Chapman-Rube sin constant 
Cj^ specific heat of solid material 

Cp constant-pressure specific heat of gas 

modulus of elasticity of SIP 


2 


f 

G 

h 

k 

L 


loc 


m 


m 


N 


Pr 


N 


Re 
p 

^base 

Pff 

^loc 

Ppeak 


modulus of elasticity of filler bar 
force on tile in x-direction 
force on tile in y~direction 
friction factor 
conductance 

unit conversion for acceleration 
tile step height 

convective heat transfer coefficient 
permeability 

lateral spring constant for SIP 
thermal conductivity 
length of duct 
tile thickness 

SIP thickness after tile installation 
filler-bar thickness after tile installation 
local Mach number in boundary layer 
local Mach number outside boundary layer 
free-stream Mach number 
mass flow rate 

mass flow rate when flow is choked 

Prandtl number 

Reynolds number 

Stanton number 

static pressure 

pressure at base of rearward-facing step 

pressure at forward-facing step over tile-to-tile gap 

local static pressure outside boundary layer 

peak static pressure at forward-facing step 
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^plateau 

^scale 


^sep 

^trough 

Poo 

P1 

P2 

R 

r 

As 

T 

T 

T 

average 


'^inlet 


loc 


T 

a 

T. 


t 

U 

u 


loc 


plateau pressure, pressure at constant pressure region 

pressure scaling factor to account for decay of plateau pressure with 
step height 

static pressure at separation point 

slight pressure drop which occurs ahead of reattachment 
free- stream static pressure 

average pressure across upstream gap of stepped tile (boundary condition 
for gap flow calculation) 

average pressure across downstream gap of stepped tile (boundary condition 
for gap flow calculation) 

gas constant 

recovery factor 

straight-line length from node to node 
static temperature 

static temperature in boundary layer 
average static temperature 

temperature of gas entering tile-to-tile gap at surface 
local static temperature 
stagnation temperature 

local stagnation temperature in boundary layer 

wall temperature 

free-stream temperature 

upstream temperature 

downstream temperature 

time 

local velocity outside boundary layer 
velocity 

local velocity in turbulent boundary layer 
average velocity 
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’^ref 

^sep 

Ax 

^ff 


z 

z 

a 

r 

A 

6 

^eff 

6* 

T) 

0 

A 

^loc 


f iller-bar-to-tile gap width 
minimum gap between filler bar and tile 
initial gap width 
tile-to-tile gap width 

Shuttle orbiter coordinates (fig. 1(b)) 

distance from leading edge 

. length of separation when h = 6 
Pf) eff 

normalized separation length 
length of separation 

longitudinal coordinate measured in downstream direction from point of 
separation 

longitudinal coordinate measured in upstream direction from a forward- facing 
step 

homogeneous solutions to energy equation for a compressible laminar boundary 
layer 

compressibility factor 

distance from wall surface 

angle of attack 

ratio of specific heats 

thermal boundary- layer thickness 

velocity boundary- layer thickness 

effective boundary-layer thickness 

displacement thickness used when flow is turbulent 

independent variable 

particle shape factor for flow through porous media 

sweep angle 

viscosity 

local viscosity outside boundary layer 
viscosity at wall 
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cor 


"pk 


m 


4 ^ 


correlated pressure ratio over extent of separation when step height equals 
boundary- layer thickness 

plateau pressure ratio over separated-f low region 

peak or reattachment pressure ratio for a given step height and boundary- 
layer thickness 

density of gas 

density of solid material 

shear at wall 

flow turning angle equivalent to a wedge angle which will produce a pressure 
rise equal to plateau pressure rise at local Mach number 


stream function 

A prime indicates a derivative with respect to an independent variable. 
Abbreviation s: 

B .L . boundary layer 

SIP strain isolator pad 

Mass flow rates: 


m 


TH 


mass flow rate per unit area in x-direction for THERM 


m 


m 


m 


TH 


FL 


FL 


ma 


ss flow rate per unit area in z-direction for THERM 


ma 


ma 


ss flow rate per unit area in x-direction from FLOW 


ss flow rate per unit area in z-direction from FLOW 


CALCULATION OF PRESSURE DISTURBANCES DUE TO STEPPED TILES 

Two locations on the lower surface of the orbiter, one on the fuselage and one 
on the left wing (see fig. 1(b) )/ were selected for analysis, since damaged filler 
bars were observed in the vicinity and flight pressure and temperature data were 
available from instrumentation at these locations. The period of interest from the 
first-flight entry trajectory (ref. 1) of the orbiter is illustrated by the curves 


shown in figure 2. The time t equals 0 when the Shuttle orbiter reaches an alti- 
tude of 400 000 ft during atmospheric entry. Prior to t = 500 sec, gap flow does 
not contribute significantly in heating the filler bar. Later than 1400 sec, flow 
velocities have decreased (M < 4) , so that the available aerodynamic heating is 
insufficient to damage the filler bar. The plotted symbols in figure 2 represent the 
trajectory times chosen for calculation of the pressure disturbances caused by a 
stepped tile. Pressure disturbances for other times in the trajectory were obtained 
by interpolation. 

The pressure and other thermodynamic properties at the edge of the boundary 
layer were calculated using a tangent-cone approximation, where the half-angle of the 
cone is equal to the local flow-deflection angle (i.e., the local body-deflection 
angle plus the angle of attack). The real-gas, axisymmetric, flow-field solution 
over a cone was obtained using a time- asymptotic numerical procedure with equilibrium 
thermodynamic properties obtained from reference 2. The tangent-cone approximation 
has been shown to yield reasonably accurate predictions of the local flow on the 
windward surface of Shuttle-like configurations. (See ref. 3.) The boundary-layer 
thicknesses used to estimate the pressure perturbations at the tile-step junctures 
were calculated by the method described in reference 4. 

The principal inviscid and viscous local flow parameters used in this study are 
presented in figure 3. Even though the flight Mach numbers ranged from approximately 
26 down to 5, the corresponding local Mach numbers for both locations on the lower 
surface of the orbiter were confined to a relatively narrow band from 3.75 down 
to 2.5. The variations in local static pressure are more dramatic, exhibiting a 
rapid rise beginning at about 700 sec and peaking near 1200 sec- 

The local boundary-layer thicknesses of interest are also presented in figure 3 
for the wing and fuselage locations. When the boundary layer is laminar, the 
boundary-layer thickness is the correlating parameter in the subsequent separated- 
flow pressure distributions. When the boundary layer is turbulent, the correlating 
parameter is the displacement thickness. The entry time at which boundary- layer 
transition from laminar to turbulent occurred was based on flight thermocouple data 
and was determined to be approximately 1200 sec for the fuselage location and 
approximately 1100 sec for the wing location. The corresponding boundary-layer 
thicknesses vary from approximately 0.87 to 4.0 in. on the fuselage and from 0.54 
to 1.54 in. on the wing. The measured postflight step heights ranged from -0.099 
to 0.12 in. in a fairly random manner; however, in this study the maximum step 
heights were limited to ±0.10 in. at the fuselage location and ±0.06 in. at the wing 
location. The ratios of step height to boundary-layer thickness were always less 
than 0.2. 


The pressure disturbances around a vertically stepped tile are illustrated in 
figure 4. The disturbance aft of the downstream face is a reduced pressure region 
where the flow leaves the surface of the stepped tile and reattaches downstream. The 
disturbance forward of the upstream face is caused by the separated flow region in 
front of the step. Both disturbances are directly over gaps between adjacent tiles. 

The base pressure on the downstream side of a vertically displaced tile is 
fairly insensitive to the range of local conditions encountered on the lower surface 
of the orbiter during entry. A base-pressure correlation for two-dimensional laminar 
and turbulent flow (ref. 5) employed in this study is presented in figure 5. The 
correlation is based on the local unit Reynolds number, — which is 

^loc 
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then multiplied by the characteristic length h 


2n-2 


. The characteristic 

length is a function of the boundary~layer thickness 6 and the step height h« The 
exponent n was empirically determined in reference 5 to be 0.9. The base-pressure 
histories of the fuselage and wing points for the STS-1 entry as determined from the 
above correlation are presented in figure 6. The two histories are almost identical. 

piT6ssure distributions ahead of the vertically displaced tile are charac— 
tBristic of a forward— facing step. If the step height is large enough/ the pressure 
distribution assumes a characteristic shape (fig. 7). The key features include the 
undisturbed local pressure a J^ise associated with boundary-layer lift-off from 

the surface, Pseo' ^ constant pressure region Ppiateau' ^ slight pressure drop 

p , , which generally occurs just ahead of the final pressure rise at reattach- 

ment. The plateau pressure also tends to be a reasonable average of the pressure 
over the extent of separation. 

The existence of a plateau pressure for both laminar and turbulent separated 
boundary layers has been experimentally verified over a wide range of flow conditions 
and for several triggering mechanisms such as forward-facing steps, flaps deflected 
into the stream, shock impingement, and compression surfaces. When the local flow 
conditions are known, the plateau pressure can be predicted fairly accurately with a 
simple algebraic equation. Several such equations have been empirically developed 
(e.g. , see refs. 6 and 7) and employed with various degrees of success. Figure 8 is 
a plot of P 1 4 . /Pi versus Mach number from test data reported in reference 8 
and e 3 q>erimental^results reported in references 6, 7, and 9. The data are directly 
applicable to the lower surface of the orbiter, since the local Mach number behind 
the bow shock ranges from 2.5 to 3.75 during the high aerodynamic heating period of 
entry. The lack of an expression which includes the influence of Reynolds number 
would normally cause doubt about the range of applicability and overall accuracy. 
Therefore, a large amount of laminar and turbulent data, including the summaries from 
two other eit^irical correlation studies, were used to validate the expression 
(refs. 6 and 7). The data in the figure fall within a ±5-percent accuracy band. 

The extent of separation of the flow ahead of a forward— facing step is dependent 
upon the local flow conditions and the step height. In this study, it was assumed 
that the dividing streamline over the extent of separation was a straight line, as 
illustrated in figure 9. The average or plateau pressure ratio over the separate d- 
f low region ^ is known, and the expression for the flow turning angle (or equiva- 
lent wedge ang?e) ^ corresponding to the plateau (or average) pressure rise can be 
found from the relationships for an oblique shock over a wedge (ref. 10). To account 
for real-gas effects, an effective ratio of specific heats y was used in the ideal- 
gas oblique shock relations at high free-stream Mach numbers. In addition, it is 
known that the boundary-layer reattachment point always occurs at the top of the 
forward— facing step. The range of validity of this assumption is dependent upon the 
nature of the boundary layer ( laminar or turbulent) and its thickness relative to the 
step height. When the boundary layer is l^inar, the step height must be equal to or 
greater than the boundary-layer thickness 5, and when the boundary layer is turbu- 
lent, the relation for x^^ shown in figure 9 is accurate only if the step height 
is equal to or greater than^1.5 times the displacement thickness 6*. These two 
parameters are designated in this study. All the step heights of concern in 

this study were less than 6gff. Therefore, an eii 5 >irical correlation of extent of 
separation for laminar and turbulent boundary layers for step heights less than 

6 was established (see fig. 10) from the data in references 11 and 12. 

ef f 
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Several experimental studies have been made of the pressure distributions ahead 
of steps whose heights are fractions of the boundary-layer thickness. (See refs. 11 
through 13.) Figures 11(a) through 11(c) present turbulent data for the range of 
local Mach numbers of interest in this study. It can be seen from these figures that 
both the plateau pressure and the character of the pressure distribution begin to 
change significantly for h/6* < 1.5, the differences becoming greater with decreas- 
ing ratios of step height to boundary- layer thickness. T!he local Mach number for the 
data of figure 11(d) is lower than any critical Mach number of this study. It is 
presented here to illustrate the effect of the limiting condition of very small step 
heights to boundary-layer thickness ratios. For a given boundary-layer thickness and 
a set of local flow conditions, the pressure distributions associated with a range of 
step heights tended to form a closely related set of curves when plotted relative to 
the separation point rather than the step face. (See ref. 13.) The information 
seems applicable to the empirical approach employed in this study, except that the 
present data set covers a wide range of boundary-layer thickness and local flow 
conditions in addition to the variation in step height. Therefore, in this study 
the experimental pressure data were normalized by calculated plateau pressures, 
and the experimental longitudinal separation lengths were normalized by the 
empirical relations for estimating extent of separation already established in 
an attempt to correlate the data. Although the data did not correlate exactly, 
an average pressure distribution for h = was identified, as well as func- 

tions for perturbing that average to represent the pressure distributions for 
step heights less than (See fig. 12.) The abscissa is simply the 

normalized separation length relative to the location of the initial pressure rise 
(i.e., X ^ \ where x., . . is calculated for h = 6 

eff eff 

and x^^^ =0 is coincident with the initial pressure rise). The curve labeled 
^cor plateau pressure ratio representing h = the pressure distribution 

which is perturbed to generate the pressure distributions at lesser step heights. 

The curve Pg^ale ^^P^^esents the decay in plateau pressure ratio with decreasing 
step height and is also used to scale the basic pressure ratio distribution ^ 
in the pressure direction. The curve ^ represents the decay in maximiim (or 
reattachment) pressure ratio with step height. 

The validity of the approach is demonstrated through the unswept forward-facing 
step data displayed in figure 11. The data presented in figure 13 are replotted from 
figure 11 relative to the separation point corresponding to a step height equal to 
^eff h/6* = 1.5 in this case). The separation point is at Ax/6* =0. A func- 

tional relationship between the family of ejgjerimental curves is quite evident and 
consistent over the range of local Mach numbers of interest. When the correlations 
developed in this study are applied to the data, the major discrepancy seems to be a 
slight underprediction of the extent of separation which is amplified at lower Mach 
numbers. However, the basic pressure distribution and the envelope of peak pressures 
seem adequate except at the lowest Mach number, at which both peak and distributed 
pressures are underpredicted. 

The data available to establish correlations for the effects of sweep on 
forward-facing steps were extremely limited. Reference 14 is a fairly comprehensive 
study of sweep effects for two step heights at two local Mach numbers (1.61 and 2.01) 
on turbulent boundary-layer separation. Even though there was an indication of a 
slight effect of Mach number on the results due to sweep, the two available Mach 
numbers were not adequate for a correlation study. The data for the higher Mach 
number were used because Mach 2.01 is closer to the Mach numbers of interest in this 
study (fig. 14). 
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Because of the limited parametric experimental data, an approach was developed 
based on the data for a local Mach number of 2.01 and the step size corresponding to 
6eff (i.e., h/6* = 1.5). It is assumed that the effects of reduced step heights 

are adequately represented by the previous correlations and that the effects of sweep 
would be adequately represented by simple ratios of swept-to-un swept pressure 
distributions at ^2.01 and h/6* = 1.5. The ratios for the four available 

sweep angles (A =0®, 15®, 30®, and 45®) are presented in figure 15. The ratio for 
A = 0® is, of course, equal to one and is shown only for reference. The decrease in 
pressure at A = 15® is nearly constant and is small over the extent of separation. 
The major difference (or reduction) in pressure for the higher sweep angles occurs 
nearest the face of the step. The 30-percent reduction in pressure for A = 30® is 
nearly constant over the 50 percent of the separation length closest to the face of 
the step. The reduction in pressure for the 45® -swept step is the primary condition 
of interest in this study, because the TPS tile pattern relative to the local flow 
direction was designed to be as near to a 45® sweep as possible. The percentage 
reduction in pressure at 45® sweep is quite large. It varies continuously across the 
length of separation and reaches a minimum value of approximately 25 percent at 
20 percent of the separation- length location in front of the step. 

The effects of these pressure reduction ratios on the unswept pressure distribu- 
tions are illustrated in figure 16 for the case when the step height is equal to 
5 The basic pressure distribution is shown in figure 12. The resulting 

distributions look reasonable, with the possible exception of the 45® sweep. How- 
ever, all the basic elements of the characteristic separated-f low pressure distribu- 
tion shown in figure 7 are retained in each distribution, with the 45° -sweep distri- 
bution being more representative of the character noted for large step sizes at lower 
Mach numbers. 


INCOMPRESSIBLE STEADY FLOW MODEL 

An existing incompressible steady flow model, originally developed for ascent 
loads analysis for the Shuttle thermal protection system (TPS) , was modified to 
include conpressibility effects in the gap. The incompressible steady flow model 
represents an array of nine tiles. The Martin Interactive Thermal Analysis System 
(MITAS), a general purpose finite-difference heat transfer program, was used for com- 
puter modeling of airflow through the TPS system. (See ref. 15.) The MITAS program 
can be used for internal incompressible flow because of the similarity between 
internal incompressible flow and heat transfer equations. 

The input to the MITAS program for steady-state flow is node and conductor defi- 
nition, or, in other words, an electrical analog representation of the geometric 
configuration. All modeling was done in three dimensions to obtain tile, SIP, and 
filler-bar internal pressures and tile-to-tile gap pressures from given surface pres- 
sure boundary conditions. 

A single-tile three-dimensional model was developed first. A sketch showing the 
flow paths through a single tile is shown in figure 17. The ceramic coating, which 
covers most of the tile surface, is considered impermeable in this study. The flow 
through the silicone-rubber bond line on the bottom of the tile is considered to be 
insignificant. Thus, areas where flow can enter or leave a tile are the terminator 
gap area and the 1/2-in. wide unbonded area on the bottom outer edge of the tile. 

Flow can enter or leave the SIP through the unbonded area on the bottom of the tile 
or through the gap between the tile bottom and the filler bar. The connections to 
the surface through which all internal flow enters or leaves are the gaps between the 
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tiles. The dimensions of the single-tile model are 6 in. by 6 in. , and the thick- 
nesses of the tile and SIP are variable. The nodal arrangement for the single tile, 
also called the subject tile in the nine-tile model, is shown in figures 18 
through 20. Figure 18 shows nodes 61 through 96, which are the nodes inside the tile 
at the surface. Figure 19 shows nodes 25 through 60 (inside the tile at the bond- 
line plane), nodes G1 through G24 (at the bottom of the gap between tiles), and nodes 
B201 through B224 (the boundary nodes for the internal flow, located at the surface 
over the gap between tiles). Figure 20 shows nodes 97 through 132, which are in the 
SIP. Conductors were defined to connect the nodes in three dimensions. Examples of 
nodes in the x-y plane that were connected by conductors would be 61 to 62, 61 to 67, 
25 to 26, 25 to 31, 97 to 98, and 97 to 103. Examples of nodes that were connected 
by conductors for flow in the z-direction would be 61 to 25, 62 to 26, 67 to 31, 

68 to 32, 25 to 97, 26 to 98, 31 to 103, and 32 to 104. There are two parallel flow 
paths from the bottom of the gap to the bottom of the tile: one flow path through 

the edge of the tile and one flow path under the edge of the tile between the tile 
and filler bar (fig. 17, view A). These flow paths are represented in figure 19 by 
electrical resistance symbols connecting gap nodes to tile bond— line nodes, such as 
nodes G2 and 26. Figure 19 also shows the conductors connecting surface boundary 
nodes to nodes at the bottom of the gap, such as the conductor shown connecting 
boundary node B202 to gap node G2. 

Using the single-tile model as a guide, an array of nine tiles was developed by 
adding eight more similar tiles to the single-tile model. Flow from one tile to 
another is through the filler bar and through the gaps between tiles. The overall 
pattern of tiles from which the multitile model was developed is shown in figure 21. 
The nine tiles in the model are shown numbered 1 through 9. The locations of the 
nodes in the multitile array are shown in figures 22 through 24. Figure 22 shows 
nodes inside the tile at the surface for each of the nine tiles. Figure 23 shows 
nodes at the tile bond line and at the surface and bottom of the adjacent gaps for 
each of the tiles. Figure 24 shows nodes in the SIP for each of the tiles. 

The flow conductance modeling for the nine-tile flow model is described in the 
paragraphs which follow. The general categories of conductance values to be deter- 
mined were tile material, SIP material, gap between tiles, tile terminator gap, gap 
between tile and filler bar, and gap between filler bar and SIP. 

The conductance for flow through tile and SIP was based on flow through porous 
media. The conductance for flow through porous media as used in the analysis is 
determined as follows: 


K pA 

c 

Conductance = — ( 1 ) 
|i As ^ ^ 


The flow through the tile in the terminator gap area is flow through porous media but 
due to the coarseness of the model it was necessary to use an entrance effect factor 
on this conductance to account for the local geometry. 

The conductance for flow through the gap between tiles and the gap between tile 
and filler bar resembles flow between parallel walls. The equations for mass flow 
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and conductance based on laminar flow between parallel walls as used in the analysis 
are as follows: 


m 


. pW^A , 

1 ^ t c dp 

12 \x ds 


2 

. pW^A 

^ . 1 ^ t c 

Conductance = — r- 

12 (1 As 


( 2 ) 


(3) 


where 


dp 

^ = Pressure gradient 


As = Distance from center of one node to the next 


The gap size and resulting conductance for flow through the gap between tiles 
was calculated iteratively for tiles which had significant side pressure forces on 
them. These side forces result from the surface pressure gradient- Figure 25 shows 
the tile-to-tile gaps in the global model which were varied. The calculations of 
tile deflection were only made on the subject tile (tile 9). The force on the sub- 
ject tile was obtained by integration of the gap pressures on the sides of the tile 
to get a resultant side force in the x- and y-directions. Neglecting rotation of the 
tile about the z-axis, the translational motion of the subject tile was then obtained 
from a spring representation of the SIP according to the following equations: 


Ax = F /K (4) 

s X s 

Ay = F /K (5) 

s y s 


where 


Ax 

s 


= Subject tile x-deflection 


F = Force in x-direction 

X 

Ay = Subject tile y-deflection 
s 


12 
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= Force in y-direction 
= Lateral spring constant for SIP 

The SIP lateral spring constant was 110 Ib/in- (ref. 16). The gaps which change are 
shown in figure 25. The shaded gaps close and were decreased from the initial value 
by the tile deflection. The cross-hatched gaps open and were increased from the 
initial value by the tile deflection multiplied by a factor to account for deflection 
of the adjacent tile. The maximum gap measured after the first Shuttle flight, 

0.126 in., was used to set this factor. This factor was 1.55 so that the gap width 
for the opening gaps would be 0.126 in. at the maximum tile deflection. A condition 
which takes priority is that the tile-to-tile gap must always be greater than the 
effective gap which exists when the two smooth surfaces are in contact. These gaps 
and the resulting flow conductance based on laminar flow between parallel walls were 
changed based on the calculated pressures until the calculated gap sizes converged. 

The gap size and resulting conductance for flow through the gap between tile and 
filler bar was calculated iteratively for all tiles in the global model. The size of 
this gap is determined by the relative normal deflection between the bottom surface 
of the tile and filler bar, and also by the silicone-rubber membrane effect. The 
size of this gap is determined in a manner similar to the tile-to-tile gap. The 
dimensions and deflections which influence this gap are shown in figure 26. The 
equation used to compute tile deflection is 


AZ 


1 



( 6 ) 


where 

Ap^ = Pressure difference between tile top surface and SIP 

A = Area which Ap acts on, loading SIP 
^1 ^ 

A^ = Area of SIP which is stressed 

= SIP thickness after tile installation 
= Modulus of elasticity of SIP 
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The filler-bar deflection is due to the pressure 
rubber surface membrane and is defined by 


differential across the silicone- 


AJl 


2 



( 7 ) 


where 


Ap Pressure difference between tile-to-tile gap and filler bar 
2 

A = Area of filler bar which Ap acts on 

P2 

^ = Area of filler bar which is stressed 

^2 

Z = Filler bar thickness after tile installation 
2 

E = Modulus of elasticity of filler bar 
2 

The desired f iller-bar-to-tile gap width is obtained based on the following 
equation: 




W 


( 8 ) 


where W is the initial gap width. The minimum limiting value for is defined 

as W • ? the effective gap width which exists when the tile bottom touches the 
f iller^bar (FB) . The final condition which was considered in the determination of 
the gap width was the silicone-rubber membrane effect. This effect causes a check 
valve action which restricts flow out of the SIP/FB more than it restricts flow into 
the SIP/FB. It was assumed that this membrane effect would dominate the gap-size 
determination. So for this analysis, it was adequate to assume that the gap height 
for flow out of the SIP/FB would be and the gap width for flow into the SIP/FB 

would be as calculated from the above equation. The tile-to-f iller-bar gaps and the 
resulting flow conductance based on laminar flow between parallel walls were also 
changed based on the calculated pressures until the calculated gaps converged. 

The flow through the gap between the SIP and filler bar was modeled as flow 
through a slot. The equations describing this type of flow are 




1/2 


( 9 ) 


14 



Conductance = C A 

D c\ Ap 


1/2 


( 10 ) 


The discharge coefficient for this slot was assumed to be 0.5. The flow cross- 

sectional area for tiles with significant side deflections is varied consistent with 
the side deflections calculated for the tile-to-tile gap deflections mentioned 
previously. 


STEADY FLOW MODEL WITH COMPRESSIBILITY 

During ascent, the gap flow was assiimed to be incompressible. However, during 
entry the pressure differences from forward- facing and rearward- facing steps were 
found to be large enough that compressibility effects would be important in the tile- 
to-tile gaps. Consequently, it was necessary to modify the existing flow model to 
account for this. 

The modification consisted of changing the tile-to-tile gap flow calculation 
from one of laminar incompressible flow to one of compressible flow, either laminar 
or transitional. The modeling of conpressible flow in the gap was accomplished under 
the assumption of conpressible internal flow of an ideal gas with constant cross- 
sectional area and frictional effects. This type of flow is often referred to as 
Fanno Line flow. The approach taken was to solve the appropriate equations of con- 
servation of mass, momentum, and energy to obtain an equation for mass flow rate as a 
function of upstream pressure and tenperature and downstream pressure. An equation 
of this form can be represented in the MITAS program as follows: 


m = f(p^ ,T^ ,p^) 


( 11 ) 


G = m/(p^ - p^) 


( 12 ) 


where 


m = Mass flow rate 


p^ = Upstream pressure 
P^ = Downstream pressure 
T^ = Upstream temperature 
G = Conductance 
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The equations used to represent this constant-area internal flow are 


Energy equation: 


T = T 

o 


.2 

m 


^22 
2p A C 
c p 


Momentum equation: 


( 13 ) 


-dp - \ = pV(dV) 

C 


Continuity; 


in = pA V 
^ c 

Ideal-gas law: 
p = pRT 

where T is the stagnation temperature (assumed to be constant). The difference 
between these equations and the equations for incompressible flow is that for incom- 
pressible flow dp and dV are zero and the ideal-gas law does not apply. 

For the case of laminar compressible flow, the momentum equation becomes 


-dp - 


^ dx = pV{dV) 
pA^W^ 


(17) 


This equation was solved simultaneously with the energy equation, continuity equa- 
tion, and ideal-gas law to obtain an equation for mass flow as a function of p-|, 
P 2 , T^, and T 2 * Initially, T 2 is assumed to be equal to T-j, The resulting 

equation for mass flow rate is 


m 


= [-b + (b^ 




(18) 
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where 


a = 



ln(pi/p2) 





2^2 
4p.A 
^1 C 


b = 12 ml 


/(‘ 


2 

c t 


c 


(0-5)RT^ 




An approximate form of the energy equation is used to adjust the initial estimate 
of T 2 * After differentiating equation (15) to obtain 


dT 


.2 

m 


p A C 

C p 


dp 


(19) 


and making the approximations 


dT » 

dp « - p^ 

P “ 


the following energy equation as used in the analysis is obtained: 


T„ = 




•2, X 

n* (P2 - Pi) 


^ 1 c p 


( 20 ) 
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where 


P2 = 


The mass flow is then reconputed using this improved value of T 2 - 

For the case of compressible transition flow, the momentum equation was obtained 
from a plot of friction factor found in reference 17* The range of transition flow 
was estimated to be 2500 < < 6000 for flow between parallel plates, or channel 

flow as it is often called- The Reynolds number as used here is defined to be 


N 


Re 


2pVW^ 2W^m 


( 21 ) 


An empirical equation for friction factor was used to represent the transition 
region. The form of the equation was 


f 

4 



+ b 


( 22 ) 


where 


f 



(23) 


The values of the constants a and b were obtained by imposing two conditions for 
a Reynolds number of 600 0. First, the value of f/4 must be 0.009. Second, the 
slope of the curve represented by equation (2 2) must equal the slope of the curve for 
the equation for turbulent friction factor. The equation for turbulent friction 
factor, as given in reference 18, is 


f 

4 


0.08 


N 


1/4 

Re 


(24) 
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The conditions to be met at a Reynolds number of 6000 are 


0.009 



+ b 


( 25 ) 


8aN 


Re 


-0.0 8N 


-5/4 

Re 


(26) 


The values of the constants a and b are obtained from these two conditions/ and 
the resulting empirical equation is 


Y= (-3.156 X 10"^Sn^ + 0.01014 

4 Re 


(2500 < N < 6000) 
Re 


(27) 


This expression for f/4 was then substituted into the momentum equation to obtain a 
transition momentum equation as follows: 


2 

8 




(28) 


Equation (28) was then solved simultaneously with the energy equation/ continuity 
equation/ and ideal-gas law, as was done for the laminar flow case, to obtain mass 
flow as a function of p^ , T-|, and T 2 « This equation for mass flow was solved 

in two iterations as in the laminar-flow case. The equation for mass flow in the 
transition region as described above was 


•4 ^-2 

am + bm + c = 0 


(29) 


and was solved for 


m to obtain the equation 



+ 



(30) 
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where 


-11 

a = (-3.156 X 10 


b = 


R Infp^/p^) 


2A^C 
c p 


2 2 

ln(p /pj r(p2 - P^) 4L 

£__L + + 0.002536 


2 2 

4p A C 
c p 


2 

A W 
C t 


q RT 
C 




The Reynolds number based on twice the gap width was always less than 6000^ so it was 
not necessary to include turbulent-flow equations. 

The mass flow was limited by sonic velocity as indicated by the expressions 


m 

s 




. 1/2 


(31) 


m < m 

s 


(32) 


where 


m = Mass flow rate for sonic velocity 
s 

Y = Ratio of specific heats 


TRANSIENT THERMAL MODEL 

The thermal model was developed to represent the insulation system and aluminum 
structure in the vicinity of a tile— to-tile gap. This modeling was done using the 
MITAS program. The region modeled was half the length of a tile-to-tile gap. The 
width of the model from the center of the gap into each adjacent tile could be set 
between 1.033 in. and 1.405 in. to change the size of the model and the location of 
nodes and boundaries when necessary. The width of the model was adequate, because 
the heat-transfer rates in the direction away from the gap are small at the boundary. 
Half the length of the tile-to-tile gap (3 in.) is adequate if the mass flow rate 
down the gap being modeled is much greater than any flow into the gap being modeled 
from the adjacent gap at either edge of the model. The top and side view of the 
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model are shown in figure 27, Figure 28 is an enlarged side view of the model- The 
parts of the tile system which were modeled were the tile, tile coating, SIP, alximi- 
num substrate, filler bar, filler-bar surface membrane, and the flowing gas in the 
gaps, SIP, and filler bar. The types of heat transfer considered were conduction, 
radiation, mass transfer, and convection. Conduction was considered in the tile, 
tile coating, SIP, filler bar, silicone-rubber membrane, and aluminum substrate. 
Radiation was considered on the tile surface and in the gap between tiles. Mass 
transfer was considered in the gaps, through the filler bar and SIP, and into or out 
of the tile. Convection was considered between the gas in the tile-to-tile gap and 
the side wall of the tile, between the gas in the gap and the filler-bar membrane, 
between the gas in the filler bar and the filler-bar membrane, and between the gas in 
the SIP and the bottom of the tile. The input to the thermal program from the flow 
program was the mass flow rate per unit area at points where mass-transfer heating 
was included in the thermal model. For the fuselage location, the boundary condi- 
tions for the thermal model were a specified temperature of the gas entering the gap 
and a tile surface heat flux versus time, from which surface temperature was calcu- 
lated. For the wing location, tile surface temperature versus time was used for a 
boundary condition rather than surface heat flux versus time. 

The numbering system used in the model is shown in figures 27 and 28. Shown in 
figure 28 are the node centers which are in the same plane. There are three planes 
like this shown in figure 27 as the 2000-, 3000-, and 4000-series locations, with the 
numbers of the nodes in figure 28 added to 2000, 3000, or 4000 for the three planes 
(for example, 2152, 3152, and 4152). The tile surface is node number 101. It is the 
same for all three planes, and therefore does not have the thousand series added to 
it. The input dimensions of the model are shown in figure 28 as RiO, R11, R12, R13, 
R20, R21, R22, R23, and R24. The width of tile nodes 2212, 2222, 2232, 2242, and 
2252 are R12/15, 2(R12/15), 4(R12/15), 8(R12/15), and 16(R12/15), respectively. The 
width of all of the tile nodes are related to R12 with this same progression. When 
the 2000- or 4000-series tile nodes are bounded by a tile corner, the nodes and con- 
ductors in the tile in that plane are changed to represent a boundary as shown in 
figure 27- 


Conduction and Radiation 

Conduction and radiation are represented in the model in the usual way for a 
finite-difference program. Conduction between nodes is represented as kA^Ai* The 
thermal conductivity k is a function of temperature, except in the HRSI insulation 
tiles where it is a function of terrperature and pressure. The alximinum substrate has 
1000- and 5000-series nodes Included in the model for edge continuity. These nodes 
are the same size as the corresponding 2000- and 4000-series nodes. Radiation from 
the surface of the tile is represented by a radiation conductor between the surface 
and a node at the deep-space teir^erature, -460®F. The only parameter which influ- 
ences this radiation from the surface is the emissivity, which was input as 0.85. 
Radiation in the gap is considered to be radiation in an enclosure formed by the 
sides of the tiles, the filler-bar membrane, a surface at the top of the gap which 
represents deep space, and surfaces at the edge of the gap to bound the model. Fig^ 
ure 29 shows the gap radiation model. Node 8001 is the top boundary and the 7000- 
and 8000-series nodes are the side boundaries. View factors for elements forming 
a right angle and having a common edge are calculated using the exact method 
(ref. 19). Others were calculated using an LaRC view factor program (ref. 20). 
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Mass Transfer 


Mass-transfer heating was of major significance only in the tile-to-tile gap, 
although it was considered in the filler bar, SIP, and through the bottom corner of 
the tile where the porous tile is not coated. Mass transfer in the thermal model 
depends on the mass flows which are calculated by the flow model. The mass transfer 
is modeled using the fluid flow conductor capability of MITAS.^ The conductance is 
the product of mass flow between the nodes and specific heat (mC^). The specific 
heat of the gas is considered to be a function of temperature. 


The actual values of mass flow are obtained by multiplying the mass flow per 
unit area, which is based on input, by the cross-sectional area for flow. For the 
case of flow in the gap, there is only one mass flow per unit area down the gap and 
one flow along the gap. These flows are calculated by the flow model, which must be 
used to estimate the mass flow rates for a finer grid of conductors in the thermal 
model. Figure 30 shows the method used to obtain the x- and z-components of mass 
flow per unit area at locations in the gap as needed in the thermal model. Mass- 
transfer conductors used the appropriate x- or z-corti>onents of mass flow per unit 
area directly, and convective conductors used the resultant mass flow per unit area 
at the location in the gap of the conductor being modeled. The mass flows of the two 
models are related as shown in the following equations {see fig. 30): 
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where the subscript i refers to node i in the thermal model, and the subscripts 
in and out refer to flow in and out of node i in the thermal model* 

The temperature of gas entering the gap from the side is set equal to the tem- 
perature of the gas in the gap in the center of the thermal model at the same depth 
from the surface. This assumption was made because the flow down the gap at the 
forward-facing step is much larger than the flow entering the gap from the side. The 
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gas at the center of the tile-to-tile gap in the thermal model is the 3000 series. 

The temperature of the gas entering the top of the gap is referred to as the gas 
i^l®t temperature. Fxgure 31 shows how mass transfer in the gap is modeled. The gas 
inlet temperature is calculated using the flow conditions at the surface and is 
discussed subsequently. 


Convection 


Convection conductors are represented in the model as the product of heat trans- 
fer coefficient h^ and surface area. Convection in the SIP and filler bar is 
modeled as forced convection through packed beds. For Reynolds numbers less than 50, 
the following equation (ref. 21) for the heat-transfer coefficient applies: 


" (0.91N 

3 p u/ Pr 
m^tn / 


-0.51 

Re 


)0 


(37) 


where 


N 


Re 


= 2 ^ 
ap0 


A value of 0.91 was the particle shape factor assumed for 0, which resulted from 
assuming cylindrical fibers. The parameter a is the solid-particle surface area 
per unit bed volume. For the gaps, the convection heat-transfer coefficient was 
calculated as follows assuming Reynolds analogy: 


N = f/8 

St Pr 


(38) 


where 


h 

N = — 

St pc u 
P 


The friction factor is determined as follows by assuming fully developed laminar flow 
between flat plates: 


f 



(39) 


where the Reynolds number is based on hydraulic diameter (ref. 22) . 
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Surface Temperature 


The tile surface boundary condition was either a specified heat flux or a speci- 
fied surface temperature. The specified heat— flux condition was obtained from a 
table in the program of surface heat flux versus time and surface -temperature. These 
data, used for the fuselage location, are shown in table I. The temperature of 
node 101, the tile surface, was calculated versus time for the specified heat-flux 
condition. 'flie specified surface temperature condition for the wing location was a 
table of temperature versus time in the program from which the temperature of node 
101 was obtained. Wing surface temperature versus time is shown in table II. 


Gas Inlet Tenperature 


The gas inlet temperature is obtained by considering the flow to be isentropi- 
cally compressed from the local pressure to the computed pressure on the forward face 
of the stepped tile: 
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(40) 


The gas average static temperature ^^average obtained from a mass flow weighted 
average of the boundary-layer static tenperature from the tile surface up to the tile 
step height h: 
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(41) 


In order to calculate the average static temperature, velocity and temperature 
distributions in both compressible laminar and turbulent boundary layers must be 
obtained- For the coitpressible laminar boundary layer, a technique developed by 
Chapman and Rubesin (ref. 23) was used to determine the velocity and tenperature 
distributions at specific times during entry. This technique approximates the tem- 
perature dependence of the viscosity of air with the following linear relationship: 


= C 


loc 


loc 


(42) 


The constant C is evaluated for each set of wall and free- stream temperatures under 
consideration. Using the viscosity relationship shown in equation (42) allows the 
boundary-layer momentum and energy equations to be independent of each other so that 
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they can be solved separately. By using a transformation that replaces the inde- 
pendent variable y with the stream function 4 ;, by nondimensionalizing all vari- 
ables, and by applying appropriate boundary conditions, the momentum equation for a 
compressible laminar boundary layer can be satisfied with the Blasius solution. Ihe 
energy equation can be written as a linear, second-order differential equation. The 
tile outer-surface temperature was considered to be uniform in the filler-bar analy- 
sis. The solution for the temperature distribution in the boundary layer is 
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(43) 


where r(r|) is the temperature recovery factor and is the function for 

n = 0 from the set of temperature distribution functions that result from 

the solution of the homogeneous part of the energy equation in the rr or z- 
direction. The functions obtained by integrating the equation 
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The recovery factor is 
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and is derived from a particular solution of the energy equation. The symbols p 
and C are arbitrary variables of integration. The term f" is one of the Blasius 
functions. The independent variable 17 is related to z in the form 
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Figure 32 shows the variations of the Blasius functions^ the recovery factor and 
integrated recovery factor, and the temperature distribution functions Y^(r|) with 
the independent variable t). Figure 33 shows the velocity and temperature distribu- 
tion in the laminar boundary layer on the Shuttle lower fuselage (X = 750 in. , Y = 
-75 in.) for a local Mach number of 2.83. 

Velocity and temperature distributions in the turbulent boundary layer were 
approximated by using the 1/7-power profiles: 
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where u is the local velocity and t is the stagnation temperature in the bound- 
ary layer. The thermal boundary layer^ thickness A was assumed to be equal to the 
velocity boundary layer thickness in the filler— bar analysis: 


A = 6 


(49) 


The perfect-gas relationships between the stagnation and static temperatures in the 
local stream and boundary layer were 


T = T 1 
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where 




- 1/2 

M = \i/(yRT) 


In order to define the static- temperature distribution in the boundary layer, equa- 
tions (50) and (51) were used to obtain the following equation (ref, 24): 



With these definitions of velocity and temperature distributions for laminar and 
turbulent boundary layers, the inlet gas temperature can be calculated. The density 
variation through the boundary layer is determined using 


^loc 

ZRT 


(53) 


where Z is the compressibility factor for air. The usual boundary-layer assumption 
that the static pressure in the boundary layer is the same as the local static pres- 
sure is used in the calculation. 


FLOW-TO-THERMAL MODEL INTERFACE 


Interfacing the flow and thermal models required processing the mass flow rates 
calculated by the flow model to obtain the mass flow per unit area needed for the 
thermal model. The flow model represents nine tiles, while the thermal model repre- 
sents only 3 in, of gap and part of the two tiles adjacent to the gap. Interfacing 
between the programs is accomplished using two subprograms. The first subprogram 
selects mass flow data from flow model results which apply to the given location of 
the smaller thermal model. The second subprogram processes the selected mass flow 
data to obtain mass flow per unit area for the thermal model (fig. 34). The input 
to the mass flow selection subprogram is information to identify the location of 
interest, and the output was the required mass flows for this location. The inputs 
to the mass flow processing subprogram are the selected mass flows and dimensions 
which are common to the flow and thermal models. The output of the mass flow 
processing subprogram is the mass flow per unit area for conductors in the thermal 
model which require mass flow information to represent convection or mass-transfer 
type of heat transfer. 
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The flow model represents steady flow at a point in time, while the thermal 
model represents a transient which spans a time interval. The thermal model requires 
flow information as a function of time, so two runs of the flow model were required 
for one run of the thermal model. One run of the flow model is made for the starting 
point of the time interval, and another run is made for the stopping point. Mass 
flow information is transferred to the thermal model for each of the two time points. 
The thermal model then uses the mass flows per unit area for these two points in time 
to obtain mass flows by linear interpolation at any time in the time interval and to 
obtain mass flows for later times by extrapolation. The flow run for the time point 
at the end of the time interval requires temperature information which must be esti- 
mated and later checked against the temperature results of the thermal model at the 
end of the time interval. 

Figure 35 is a flow chart of the tile flow and thermal analysis. Runs of the 
flow program are indicated by FLOW, and runs of the thermal analysis program are 
indicated by THERM. THERM uses the steady-state mass flows and other temperature 
boundary conditions to calculate the transient thermal response of the gap region. 

The input temperatures are known for the run of FLOW at the start of a time interval, 
but the input temperatures for the run of FLOW at the end of a time interval must be 
estimated and later checked against temperatures which are calculated by THERM. If 
the input temperatures estimated for the last run of FLOW are more than 5 percent 
different from the calculated temperatures, the last run of FLOW must be repeated 
with new estimated input temperatures, and the last run of THERM must also be 
repeated using the new set of mass flows calculated by FLOW. The thermal model THERM 
calculates for two time intervals per run. The results of the first time interval 
are saved, and the results of the second time interval are used to estimate input 
temperatures for the next run of FLOW. 

The desired location of the thermal model in the tile array to examine filler- 
bar charring was the location where the highest gap temperatures would occur. High 
gap temperatures are caused by large mass flow rates down the gap, and therefore the 
thermal model location in the tile array was the part of the gap which had the high- 
est calculated mass flow rate down the gap. The calculated mass flow rates down the 
gap for runs of the flow model with a stepped tile were always greatest at the same 
gap location, so this location was chosen for the thermal model. The location of the 
thermal model thus determined is shown in figure 36- The necessary condition for 
accuracy of the thermal model, that the mass flow rate down the gap be much greater 
than the mass flow rate along the gap, was satisfied at this location. 


CALCULATED FILLER-BAR TEMPERATURES DUE TO STEPPED TILES 

A thermal analysis was performed for two locations on the bottom surface of the 
Shuttle orbiter for the STS- 1 entry conditions. The two locations were a fuselage 
point (500 in. from the nose) at orbiter coordinates of X = 750, Y = -75, and 
a wing point (105 in. aft of the local leading edge) at X = 1220, Y = 275 
(fig. 1(b)). These two points were selected because (1) filler-bar charring 
occurred, (2) they are in regions where the local aerodynamic surface flow conditions 
could be computed, and (3) they represent two extremes in terms of location relative 
to the leading edge. Table III shows the extremes in step and gap dimensions of STS- 
1 postflight ground measurements made at damaged filler-bar locations. The table 
also shows definitions of filler-bar charring categories. This information was 
supplied by the Johnson Space Center (JSC). The gap heating analysis model was based 
upon the STS-1 entry trajectory for various step heights and nominal gaps, starting 
from an initial temperature of 40®F. 
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Results were obtained for tile steps of 0.1 in. and 0.04 in. at the fuselage 
location. Figure 37 shows how the ratio of step height to boundary-layer thickness 
h/6 varies at the fuselage location X = 750, Y = -75 for a constant step height 
of 0.1 in. The thickness of the tile at this location was 1.41 in., and the initial 
tile-to-tile gap size was 0.05 in. Figure 38 shows the forward and aft face pres- 
sures for the 0.1-in. stepped tile along with the local static pressure used in the 
analysis. Figure 39 shows the temperature response of the filler bar at the fuselage 
location, along with the local static gas temperature, the inlet gas temperature, and 
the tile surface temperature. Flight data indicated that boundary-layer transition 
from laminar to turbulent flow starts at 1235 sec* The inlet gas temperature 
increases at this time due to the change in boundary-layer temperature profile. A 
peak filler-bar temperature of 1502®F occurred at t = 1285 sec. Figures 40 and 41 
give a similar set of results for 0.04-in. tile step height. The calculated tempera- 
tures for the filler-bar membrane node which reaches the highest temperature (node 
4022) are given in figure 42 for the fuselage location. For this case, the filler- 
bar temperature peaks at 969°F. These two step height runs at the fuselage location 
represent category 3 and category 1 charring conditions (table III). The filler-bar 
temperature distribution at the time when the peak temperature occurs is given in 
table IV. 

Results were obtained for three tile steps at the wing location X = 1220, 

Y = 275. The tile thickness at this location is 2.25 in. The step heights which 
were run are 0.017 in. , 0.03 in. , and 0.06 in. The 0.017— in. step height case was 
run because 0.017 in. is the maximum design step height allowed at this wing loca- 
tion. The 0.03-in. and 0.06-in. step height cases were run to represent category 1 
and category 3 filler-bar charring, respectively. Figures 43 through 48 show the 
pressures and calculated temperatures for the wing location. The peak filler-bar 
temperatures calculated for the three step heights were 849«F for the 0.017-in. step, 
1284®F for the 0.03-in. step, and 1748®F for the 0.06-in. step. The calculated tem- 
peratures for filler-bar membrane node 4022 are given in figure 49 for the wing loca- 
tion. The calculated peak temperatures and temperature distribution for the filler- 
bar membrane are given in table IV. The calculated gap-width variation during entry 
for representative cases is given in table V. The results indicate that filler-bar 
temperature is strongly dependent on step height. 

The effect of early boundary-layer transition on peak filler-bar temperature was 
studied by assuming transition at 1225 sec (10 sec earlier) for the 0.1-in. step case 
at the fuselage location. The calculated filler-bar temperature peaked at 1578®F, 
compared with 1502®F for transition at 1.235 sec. 


POTENTIAL SOLUTION TO THE PROBLEM OF FILLER-BAR CHARRING 

The previous analyses have shown that filler-bar charring can be caused by the 
flow of hot gases in the tile-to-tile gaps as a result of forward-facing and 
rearward-facing steps at the surface. The flow of hot gases in the gap can be 
reduced by installation of gap fillers between the tiles with excessive steps. Many 
ceramic fiber products have the combination of high flow resistance and high 
temperature capability, which is needed for a gap filler. 

The thermal model was used to estimate the performance of a half-height gap 
filler. A sketch of the gap filler is shown in figure 50. The material chosen for 
the gap filler was Carborundum Fiberfrax 880 JH. The density of this material is 
6 Ibm/ft .The permeability was measured experimentally and found to be 
3.66 X 10 ft^ at a pressure of 142 psf (pounds per square foot). The pressure 
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gradient was 1114 psf per foot. (See table VI for test results.) This gap-filler 
material has a tenperature capability of 2600^F. The adjustment to the model to 
represent the gap filler was to set the tile-to-tile gap such that the flow resist- 
ance of the gap was equivalent to the flow resistance of a 0.05-in. gap with a half- 
height ceramic fiber gap filler installed. First/ the size gap which is equivalent 
to a 0.05-in. gap with a full gap filler installed was determined. Based on a gap- 
filler permeability of 3.66 x 10 ft^z the equivalent gap width is 0.00316 in. An 
approximation for the equivalent gap width for a half-height gap filler of this type 
is obtained by multiplying the full gap-filler equivalent gap width by the cube root 
of 2. The half-height gap-filler equivalent gap thus obtained is 0.00398 in. The 
thermal model was run with an initial gap width of 0.004 in. to represent the instal- 
lation of this half-height gap filler in the 0.05-in. tile-to-tile gap. The wing 
location with a 0.06-in. tile-to-tile step was chosen for the gap-filler run analy- 
sis. The pressure boundary conditions used for this run are the same as for the 
analysis without the gap filler and are given in figure 31. The temperature results 
of the gap-filler analysis are given in figure 51. The results show that the addi- 
tion of the gap filler protected the filler bar, so that charring would not occur. 

The maximum calculated temperature of the filler— bar membrane was only 450®F / while 
970®F is required to start the charring process. The calculated filler-bar tempera- 
ture for the same case, but with a 0.05-in. gap without a gap filler/ was 1748°F. 

The gap filler reduced the filler-bar temperature by 1298°F. 


CONCLUDING REMARKS 

An analysis of the heating that occurs during entry in the tile-to-tile gap due 
to uneven tile heights in the Shuttle orbiter thermal protection system has been 
performed. In order to conduct the analysis/ a technique was developed for predict- 
ing the local pressure disturbance caused by a stepped tile. In addition/ a method 
was developed for calculating the airflow rates in the tile— to— tile gapS/ and for 
predicting the temperature level of the gas ingested from the boundary layer. A 
thermal analysis was conducted on the stepped-tile configuration to determine the 
extent of heating on the tile sidewalls/ the filler bar/ and the SIP due to the hot- 
air flow in the gap. Combinations of tile step heights and resulting tile-to-tile 
gaps that could cause varying degrees of damage to the filler bar on the lower fuse- 
lage and wing were determined. The magnitudes of the predicted step heights and gaps 
were comparable to the step heights and gaps that were observed in damaged regions 
after the initial flights of the Shuttle. Filler-bar charring may be prevented by 
controlling steps within tight tolerances during installation and maintaining these 
tolerances in flight. The acceptable in-flight steps were 0.04 in. at the fuselage 
location and 0.017 in. at the wing location. If the tolerances cannot be maintained/ 
tile-to-tile gap filler would be an alternative. A ceramic fiber gap filler was 
investigated analytically using the techniques described in the paper. The results 
indicate that the gap filler would reduce the airflow rates sufficiently to prevent 
damage to the filler bar. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton/ VA 23665 
July 27/ 1983 
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TABLE II." SURFACE TEMPERATURE DATA USED FOR WING POINT 

X = 1220, Y = 275 

From Space Shuttle Orbiter Entry Aerodynamic Heating 
Data Book, Rockwell Publication SD73-SH— 0 184 , 
Revision C, October 1978 


Time, sec 

Wall temperature, °F 

0 

40 

505 

1830 

634 

1820 

811 

1815 

915 

1805 

1013 

1720 

1 104 

1540 

1200 

1478 

1298 

1630 

1335 

1542 

1375 

1423 





TABLE III 


STS-1 MEASUREMENTS IN DT^GED FILLER-BAR AREAS 


STS-1 PostfXight measurements in damaged filler-bar areas 
Tile-to-tile steps -0-099 in. to 0.12 in. 

Tile-to-tile gaps 0.010 in. to 0.13 in. 


Category 

1 

2 

3 


Categorization of damaged filler bar 
Temperature Damage 

970*^ F Silicone-rubber membrane 

1100®F Silicone-rubber membrane 


1375°F 


Silicone-rubber membrane 
bar charred 


discolored 
charred 
and filler 


Approximately 600 damaged filler bars on STS-1 






































TABLE VI.- PERMEABILITY TEST RESULTS 

FOR CARBORUNDUM FIBERFRAX 880JH 

CERAMIC FIBER 

BLANKET 


Sample thickness 

= 0.093 in. 



Sample width 

= 1.26 in. 



Sample length 

= 1.26 in. 



Sample frontal area 

= 0.1172 in^ 



Standard temperature = 70®F 



[standard pressure 

- 2116 psf J 



Flow, 

standard ft^/min 

Pi, 

psf 

P2^ 

psf 

m, 

Ibm/sec 

dp, 

ds' 

psf/f t 

ft^ 

3.53 X 10"3 

200.5 

83.5 

4.409 X 10“® 

1114 

3.66 X 10"''° 

7.06 

273.4 

83.5 

8.818 X 10“® 

1816 

3.57 

10.6 

338.3 

83.5 

1.323 X 10*^ 

2426 

3.40 

14. 1 

396.8 

83.5 

1.764 

2893 

3.24 

17.6 

448.3 

83.5 

2.205 

3474 

3.14 

14. 1 

398.5 

83. 5 

1.764 

2999 

3.21 

10.6 

342.5 

83.5 

1.323 

2466 

3.31 

7.06 

278.5 

83.5 

8.818 X 10”^ 

1856 

3.45 

3.53 

203.3 

83.5 

4.409 X 10"^ 

1141 

3.54 


















Time, sec 

Figure 2.- Reconstructed STS- 1 entry trajectory and approximated values used 

in filler-bar charring analysis. 
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Figure 4.- Schematic of flow patterns and key parameters associated with 

vertically displaced tile. 
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Figure 9,- Method for estimating length of separation 







Figure 11.- Experimental separated- flow pressure distributions for turbulent flow ahead of 
forward-facing steps whose heights are fractions of boundary-layer thickness. 
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Figure 11.- Concluded 








Figure 13.- Con^arison of correlation and experimental pres sure- profile data for 
separated flows ahead of forward-facing steps. 






Figure 13* - Continued, 



Figure 13,- Continued 




Correlation curve for 
envelope of peak pressures 



Figure 1 3. - Concluded, 




ahead of un swept forward- facing steps. 
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Figure 22,- Tile array showing surface nodes (flow model)- 
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Figure 23.- Tile array showing bond-line nodes (flow model). 
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Figure 25.- Nine-tile array showing gap changes due to x- and y-def lections of 

subject tile, tile 9 (flow model) . 
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Figure 26 •- TPS vertical dimensions and deflections due to pressure loads. 
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Figure 31.- Control volume nodes for mass transfer in gap. 
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Figure 33.- Velocity and temperature distributions in laminar boundary layer on 
orbiter fuselage at location X = 750 in. , Y = -75 in. 




From (t) to (t + At) 


Figure 34.- Flow chart for interface from FLOW to THERM. 
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Figure 36.- Location of thermal model boundaries- 




Figure 37,- Ratio of step height to boundary-layer thickness at fuselage location 

for 0,1-in, step. 




Local pressure 



Figure 38.- Pressures at fuselage location for 0.1 -in. step. 



Figure 39-- Gap temperatures at fuselage location for 0.1-in. step and 

0.05-in. initial gap. 




A Local pressure 
O Forward-face pressure 



Figure 40,- Pressures at fuselage location for 0.04-in. step. 
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.05-in. initial gap. 


Category 3 



Figure 42.- Calculated temperatures at fuselage location for filler-bar membran 

node of highest peak temperature, node 4022. 
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Figure 44<- Gap temperatures at wing location for 0,017-in. step and 

0,05-in. initial gap- 
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A Local stream temp 
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Figure 46.- Gap teicperatures at wing location for 0 -03-in. step and 

0. 05-in. initial gap. 
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Figure 47,- Pressures at wing location for 0,06-in- step. 
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Figure 48.- Gap temperatures at wing location for 0.06-in. step and 

0.05-in. initial gap. 
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Figure 49*- Calculated temperatures at wing location for filler-bar membrane 
node of highest peak temperature, node 4022- 
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